Bose-Einstein statistics in thermalization and photoluminescence of quantum well 
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Quasi-equilibrium relaxational thermodynamics is developed to understand LA-phonon-assisted 
thermalization of Bose-Einstein distributed excitons in quantum wells. We study the quantum- 
statistical effects in the relaxational dynamics of the effective temperature of excitons T = T(t). 
When T is less than the degeneracy temperature To, well-developed Bose-Einstein statistics of 
quantum well excitons leads to nonexponential and density-dependent thermalization. At low bath 
temperatures Tt — > the thermalization of quantum-statistically degenerate excitons effectively 
slows down and T(t) oc l/\nt. We also analyze the optical decay of Bose-Einstein distributed exci- 
tons in perfect quantum wells and show how nonclassical statistics influences the effective lifetime 
T op t- In particular, r op t of a strongly degenerate gas of excitons is given by 2tr, where tr is the 
intrinsic radiative lifetime of quasi-two-dimensional excitons. Kinetics of resonant photolumines- 
cence of quantum well excitons during their thermalization is studied within the thermodynamic 
approach and taking into account Bose-Einstein statistics. We find density-dependent photolumines- 
cence dynamics of statistically degenerate excitons. Numerical modeling of the thermalization and 
photoluminescence kinetics of quasi- two-dimensional excitons are given for GaAs/AlGaAs quantum 
wells. 
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I. INTRODUCTION 

The fundamental features of relaxation and photolumi- 
nescence (PL) of the excitons in quantum wells (QWs) 
originate from the quasi-two-dimensionality (quasi-2D) 
of the system. In the last decade these basic processes in 
GaAs QWs have attracted continual attentionEfO. The 
formation of QW excitons through LO-phonon cascade 
emission, LA-phonon and carrier-carrier scattering com- 
pletes withjjiJZO ps after the initial excitation of electron- 
hole pairsofj'H The created hot QW excitons then ther- 
malize through low-energy acoustic phonon scattering. 
This process occurs in a sub-ns time scale. In GaAs QWs 
the relaxation of hot excitons can be observed through 
resonant PL from the optically-active bulk modes which 
refer to a small in-plane momentum. Thus the rise and 
decay times of excitonic PL relate to the thermalization 
process. Furthermore, the PL kinetics of long-lifetime in- 
direct excitons in a high-quality GaAs/AlGaAs coupled 
quantum well (CQW) is now the subject of experimental 
studyEia. There has been recent progress in the inves- 
tigation of thermalization of quasi-2D excitons in ZnSe 
QWslii In the latter case the picosecond LA-phonon- 
assisted kinetics of QW excitons is visualized through 
LO-phonon- assisted PL for all in-plane modes pn. 

While the importance of interface disorder and the 
localization affects were recognized in the very first 
experimental and theoretically studies on the relax- 



ation kinetics of QW excitons, the quality of GaAs QWs 
continuously improves from^n inhomogeneous excitonic 
linewidth of about 7.5 r^VEl toward a homogeneous one 
between 80 and 120 /xeVlj. In the present work we inves- 
tigate LA-phonon-assisted relaxation kinetics in a per- 
fect QW emphasizing the importance of Bose-Einstein 
(BE) statistics of quasi-2D excitons. In particular, we at- 
tribute jthe density-dependent PL kinetics reported, e.g., 
in RefsJ-3'0, to nonclassical statistics of QW excitons. In 
the past, most theoretical, modeling of the relaxation 
kinetics in QWsafflli3'L§li3 dealt with a classical gas 
of Maxwell-Boltzmann distributed excitons. Quantum- 
statistical effects were included in numerical simulations 
of the relaxation kinetics of a trapped quasi-2D exciton 
gasliZl and in the study of exciton-biexciton law of mass 
action in QWsEj. 

Crossover from classical to quantum statistics occurs 
near the degeneracy temperature Xn, given by 




P2D 



(1) 



where M x is the in-plane translational mass of a QW 
exciton, p2D is the 2D concentration of excitons, and 
g is the spin degeneracy factor of QW excitons. We 
put g = 4, because fox. GaAs QWs the exchange inter- 
action is rather weakEj. Furthermore, this interaction 
is completely suppressed for indirect excitons in CQWs. 
One has a classical gas of QW excitons at temperatures 
T 3> To. BE statistics smoothly develops with decreas- 
ing T ~ T leading to occupation numbers AT p > 1 of 
the low-energy in-plane modes py. We will show that 
at helium temperatures nonclassical statistics affects the 
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thermalization process already at moderate densities of 
QW excitons p 2 D>3-5x 10 10 cm" 2 . 

In this paper we develop relaxational thermodynam- 
ics and study within this approach how BE statistics in- 
fluences thermalization and photoluminescence of quasi- 
2D excitons. Relaxational thermodynamics requires 
that exciton-exciton interaction is much stronger than 
exciton-LA-phonon coupling and is appropriate if the 
concentration p2D of QW excitons is larger than some 
critical density p^d- In this case, the the exciton system 
establishes a quasi-cquilibrium temperature T. For GaAs 
QWs we will get an estimate p\ D ~ 1 — 3 x 10 9 cm" 2 . 
Equation (pi]), which is the basic equation of the re- 
laxational thermodynamics of QW excitons, provides us 
with an unified description of the thermalization process. 
The thermodynamic approach yields temporal evolution 
of the effective temperature of QW excitons T = T(t) 
from the initial value Tj = T(t — 0) to the bath temper- 
ature T(,. While we study the thermalization process in 
the both limits, classical and degenerate, of a gas of QW 
excitons, the most interesting results refer to T < T . In 
this case one finds a density-dependent nonexponential 
relaxation of quasi-2D excitons. Both acceleration and 
slowing down of the thermalization kinetics may occur 
due to BE statistics. However, with the decreasing bath 
temperature X& slowing down of the relaxation starts to 
dominate in a quantum gas. In particular, for Xj, — > we 
derive 1/lnt cooling law for QW excitons. 

The thermalization kinetics of excitons in GaAs QWs 
can be observed through resonant PL from the in- 
plane cadiative modes. Therefore we generalize the PL 
theorymij to well-developed BE statistics of QW excitons 
and model numerically the PL process in GaAs QWs. 
It will be shown that at low temperatures nonclassical 
statistics changes the law T opt oc T, where r opt is the ef- 
fective radiative lifetime of QW excitons. Furthermore, 
we calculate the temperature-independent corrections to 
the classical linear behavior of T opt (T). The corrections 
originate from BE statistics and can be traced even at 
high temperatures. 

Relaxational thermodynamics refers to the phonon- 
assisted thermalization kinetics of QW excitons. We will 
analyze the relaxation kinetics due to bulk LA-phonons, 
assuming the initial distribution of hot QW excitons is 
below the threshold for optical phonon emission. The 
bulk LA-phonons are due to a semiconductor substrate 
of a QW. While the lifetime of long wavelength acous- 
tic phonons are in a sub-/is time scale, the scattering 
LA-phonons leave and enter the QW area within a time 
~ L z /v s ~ 1 — lOps (L z is the thickness of a QW, v s is the 
longitudinal sound velocity). Therefore, we assume the 
Planck distribution of the LA-phonons interacting with 
QW excitons. The LA-phonon-assisted kinetics will be 
treated for an isolated band of ground-state QW excitons. 
Being applied to heavy-hole excitons in GaAs QWs with 
L z < 100 A this assumption means QW exciton energies 
E < 10 meV and temperatures less than 100 K. Only 
the exciton-LA-phonon deformation potential interaction 



is included in our model. 

The in-plane momentum is conserved in the scattering 
of QW excitons by bulk LA-phonons. The momentum 
conservation in the z-direction (the QW growth direc- 
tion) is relaxed. As a result, a scattered QW exciton 
interacts with a continuum spectral band of scattering 
LA-phonons of a given direction of propagation. In con- 
trast, an exciton in bulk semiconductors couples in Stokes 
or anti-Stokes scattering only with one phonon mode of 
a given direction. In Fig. 1 we depict a schematic pic- 
ture of the states which are involved in the scattering of 
a QW exciton with in-plane momentum p|| = 0. The en- 
ergy state E = couples to the continuum energy states 
E > E = 2M x v 2 , i.e., to the QW states which lie in- 
side the acoustic cone given by E = E(pu) = hv s p\\. 
The energy Eq is an important parameter of the re- 
laxational thermodynamics of QW excitons. For GaAs 
QWs with M x — 0.3 mo (mo is the free electron mass) 
and v s — 3.7 x 10 5 cm/s one has E = 46.7 peV and 
E /k B = 0.54 K. 
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FIG. 1. Schematic picture of the energy-momentum con- 
servation in scattering of QW excitons by bulk LA-phonons. 
The parabolic and conical surfaces refer to the excitonic 
E — h 2 p^/2M x and acoustic E = hv s p\\ dispersions, respec- 
tively. The bold rings show the states hp\\ > 2M x v 3 of QW 
excitons which couple to the ground-state mode py = 0. 

The relaxation kinetics of QW excitons coupled to 
thermal bulk LA-phonons is given by the Boltzmann 
equation (see, e.g.j23): 
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where the dimensionless energy is s = E/Eq = 
E/(2M x v 2 ), nf = l/[exp{eE /k B T b ) - 1] and N e are 
the distribution functions of bulk LA-phonons and QW 
excitons, respectively, is the Heaviside step-function. 
The integration variable u is given by u — cos 9, where 
< 6 < 7r is the angle between z-axis and the wavevec- 
tor of a scattering bulk LA-phonon. The scattering 
time is t sc — (n 2 h i p)/(D 2 M x v s ), where p is the crys- 
tal density and D is the deformation potential. The 
form-factor F z (x) = [sin(x)/x][e zx /(l t— j^ 2 / 71 " 2 )] refers to 
an infinite QW confinement potentiaEH. This function 
describes the relaxation of the momentum conservation 
in the z-direction and characterizes a spectral band of 
LA-phonons, which effectively interact with a QW exci- 
ton. The dimensionless parameter a ~ 1 is defined by 
a = (L z M x v s )/h. 

The kinetic Eq. (|2|) deals with an isotropic in-plane 
distribution of QW excitons, i.e., the occupation number 
iVp,| of the in-plane mode p|| relates to the distribution 
function Ne by AT P|| = Ne p = Ne- This approximation 

corresponds to the experimental conditions of Refs. [1- 
10]. The first (second) term in the figure brackets on the 
right-hand side (r.h.s.) of Eq. (||) describes the phonon- 
assisted Stokes (anti-Stokes) scattering. In accordance 
with Eq. (||), the relaxation kinetics into the ground-state 
mode p|| = is given by 
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The integral on the r.h.s. of Eq. (Q) characterizes the 
coupling of the ground-state mode E = to continuum 
of the energy states E > E (see Fig. 1). 

In Sec. II, we develop the relaxational thermodynam- 
ics of QW excitons coupled to bulk LA-phonons. A ba- 
sic equation for the effective temperature T of quasi- 
equilibrium QW excitons is derived from LA-phonon- 
assisted kinetics given by Eq. (0). The conditions of the 
validity of the thermodynamic picture are analyzed and 
tested for excitons in GaAs QWs. 

In Sec. Ill, the thermalization law T = T(t) is stud- 
ied for both classical (T ^> T ) and well-developed BE 
statistics (T <C To) of quasi-equilibrium QW excitons. 
We demonstrate that BE statistics strongly influences 
the thermalization process and leads to the density- 
dependent characteristic thermalization time r t h — 
T~th(p2D) for \T — Tf,\ <C Tb and to the nonexponential 



density-dependent relaxation at T — Tb > Tb. In particu- 
lar, for low bath temperatures Tb < E$/kB we find a very 
slow thermalization T(t) oc 1/ In t of quantum degenerate 
quasi-2D excitons. The numerical simulations of the re- 
laxational dynamics are given for excitons in GaAs QWs. 

In Sec. IV, we develop a theory of resonant PL of sta- 
tistically degenerate QW excitons. An effective radia- 
tive lifetime r opt of a quasi-2D excitonic gas is calculated 
for BE-distributed QW excitons. We show that the law 
T opt cx T, which is valid for classical statistics of QW ex- 
citons, is violated at low temperatures and at T — ► the 
effective decay time T opt tends to 2rjj, where tr is the 
intrinsic radiative lifetime of a QW exciton. The PL ki- 
netics of QW excitons is described in the thermodynamic 
approach within three coupled equations for T(t), p2£>(t), 
and r op t(T, Pid)- We give numerical modeling of the T- 
and p2D-dependent PL kinetics of quantum degenerate 
quasi-2D excitons. 

In Sec. V, we discuss the influence of interface polari- 
tons and QW biexcitons on the relaxation kinetics of sta- 
tistically degenerate excitons in perfect QWs. Our the- 
ory is always appropriate for thermalization of indirect 
excitons in CQWs, where even at small densities p2D 
the interface polariton effect is rather weak and biexci- 
ton states are unbound. For single QWs we argue that 
QW polaritons and biexcitons are considerably weakened 
at large densities, due to particle-particle scattering. In 
this case our model is applicable to usual QWs with di- 
rect excitons. 

In Appendix A, in order to estimate the critical density 
p\ D of QW excitons for the development of relaxational 
thermodynamics we calculate the characteristic equili- 
bration time t x _ x due to exciton-exciton scattering. The 
equilibration rate 1/t x - x is found in both limits of clas- 
sical and well-developed BE statistics. 

In Appendix B, we analyze how a relatively large ho- 
mogeneous linewidth ~ Ti/t x - x 3> h/r^h, ^/tr of high- 
density QW excitons in a perfect QW influences the PL 
process. 



II. RELAXATIONAL THERMODYNAMICS OF 
QW EXCITONS 

In this section we summarize the thermodynamic re- 
lations for an ideal two-dimensional gas of bosons and 
derive the basic equation for relaxational dynamics of 
QW excitons. 



A. Thermodynamic relations for quasi-2D excitons 

The thermodynamic equation p = p{T,p2D) for a 
quasi-equilibrium gas of QW excitons can be derived 
from the condition: 
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From Eq. M) one obtains 



p = k B T In (l-e" To/T 



(4) 



(5) 



where the degeneracy temperature To is given by Eq. (|l]) 
with (7 = 4. With the chemical potential pL of Eq. (||) the 
equilibrium distribution function of QW excitons is 



JV S e E/k B T + e -T /T _ 1 



(6) 



In particular, the occupation number of the ground-state 
mode is given by N^ =0 = N e E q =0 = exp(T /T) - 1. 

Classical, Maxwell-Boltzmann, statistics of QW exci- 
tons is realized for T » T . In this case Eqs. 
reduce to 



p = fcsTln(To/T) and 



N eq 



T /T « 1 , (7) 



i.e., the occupation numbers of QW modes p|| are much 
less than unity. In the opposite limit T < Tq of well- 
developed BE statistics one has 



-k R Te- To/T and 



N eq 



o T /T 



> 1 



(8) 



According to Eq. (0), the chemical potential of a degen- 
erate gas of QW excitons approaches zero much faster 
than the temperature T. While the chemical potential 
/i < for T > and the BE condensation of the QW 
excitons is absent within the thermodynamic approach 
for temperatures above zero, the occupation number of 
the ground-state mode p|| increases exponentially with 
decreasing T <C T . 

Equations (§)-(|) can be applied to an arbitrary 
two-dimensional quasi-ideal gas of Bose-particles with 
quadratic dispersion. The specific characteristics of the 
bosons, like the spin degeneracy factor g and the transla- 
tional mass M x , enter the thermodynamic relationships 
only through the degeneracy temperature To defined by 
Eq. ®. 

For GaAs QWs one estimates T = (ttH 2 p 2 o) / (2M x k B ) 
= 4.6 K {k B T = 399.5 peV) for p 2D = 10 11 cm" 2 . This 
density of quasi-2D excitons corresponds to the mean in- 
terparticle distance ~ 1/ ^/p 2 D — 320 A and to the Mott 

parameter P2d[o^x D " 1 ] 2 — 0.04, i.e., the QW excitons are 
still well-defined quasiparticles. Here a x 2D ^ is the Bohr 
radius of a quasi-2D exciton. In next subsection we will 
use Eqs. (§|)-(j|) to develop the relaxational dynamics. 



B. Thermalization equation for quasi-equilibrium 
excitons in QWs 

The thermalization of QW excitons occurs through 
a sequence of quasi-equilibrium thermodynamic states, 



which are characterized by an effective temperature T = 
T(t) and chemical potential p — p(t), provided that the 
exciton-exciton interaction in a QW is much stronger 
than the coupling of QW excitons with bulk LA-phonons. 
The initial interaction is conservative and equilibrates the 
system of QW excitons without change of total energy. 
The characteristic equilibration time t x _ x depends on 
the density p 2B of QW excitons. For excitons distributed 
below the threshold for LO-phonon emission, energy (ef- 
fective temperature) relaxation results from QW exciton 
- bulk LA-phonon scattering and is characterized by an 
effective thermalization time r t h- The hierarchy of in- 
teractions means that in a large and interesting range of 
density we have t x - x -C r t h- In this case the equilibra- 
tion and thermalization kinetics can be naturally sepa- 
rated (see, e.g., Ref.E3) and the thermodynamic approach 
is correct. 

The condition t x - x -c fth will hold for p 2 u > P 2 d, 
where p 2D is some critical density for QW excitons. 
For a classical distribution of quasi-2D excitons one has 
T /T < 1 and N eq = (T /T) exp(-E/k B T) < 1, ac- 
cording to Eq. (0). In this limit the characteristic equi- 
libration time t x ^ x is estimated by Eq. (A4) of Ap- 
pendix A, while the thermalization time r t h is given by 
Eq. (15) of section III. The comparison of Eq. (A4) and 
Eq. (15) yields t x - x <^ Tth provided that 



4h I p, x 



a 



2D 



(9) 



where p x is the reduced mass of a QW exciton and the 
constant C 2 u — C 2 o(a) ^ 1 is defined by Eq. (|l5|). 
For a classical gas of excitons in GaAs QW of thickness 
L z — 100 A we estimate p 2D = 1.2 x 10 9 cm~ 2 and the 
corresponding temperature scale Tq = 56 mK . For com- 
parison, the critical density of|-about 4 x 10 9 cm~ 2 has 
been estimated in experimentsu with high-quality GaAs 
QWs of L z = 80 °A. For well-developed BE statistics of 
QW excitons one has To /T 3> 1 and the occupation num- 
bers of the energy states E > k B T exp(— Tq/T) are given 
by the Planck function, i.e., N E q >0 = l/[exp(E/k B T)-l], 
according to Eq. (^|). In this case t x - x and r t h are given 
by Eq. (A5) of Appendix and Eq. ( |lS| ) of section III, re- 
spectively The condition t x _ x <ti r t h reduces to 
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T 



3 T /T 



4/v 



(1 - 7r/4) IM, 




2D 



(10) 



where the constant C2D — C2d(o) S> 1 is given by 
Eq. ([l8]). The inequality ( ]To| ) always holds for strongly 
degenerate QW excitons. 

The criteria (^) and ( ]l0| ) of the validity of the ther- 
modynamic picture at T > To and T <C T), respec- 
tively, are independent of the scattering length ~ a x 2D ^ 
of exciton-exciton interaction. This is due to the quasi- 
two-dimensionality of QW excitons (for details see Ap- 
pendix A). In contrast, in a three-dimensional gas the 
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equilibration time t x - x due to particle-particle, interac- 
tion explicitly depends on the scattering lengthen. In fur- 
ther analysis we assume the hierarchy of interactions, i.e., 
that p2D > P2D so that the inequality (|J) is valid. The 
initial density of photogenerated excitons in GaAs QWs 
is usually larger than 10 9 cm" 2 [3-10], indicating that 
the thermodynamic picture of relaxation is adequate for 
typical experimental conditions. We will omit the super- 
script in N E q , because within our approach Ne = N E q . 

The thermalization dynamics of QW excitons is de- 
termined by the slowest elementary LA-phonon-assisted 
relaxation process of the kinetic Eqs. (||) and (||). The 
joint density of states for Stokes and anti-Stokes scatter- 
ing from the energy mode E continuously decreases with 
decreasing E, according to the r.h.s. of Eq. (||) (note, 
that the density of quasi-2D excitonic states is constant 
given by 16n 2 M x /Ti 2 for g = 4). Moreover, the low- 
energy states E < Eq/ '4 are not active in Stokes scat- 
tering and couple with the corresponding LA-phonon- 
separated states E > E only through the anti-Stokes 
process. The above arguments show a particular status 
of Eq. (^), which describes the relaxation kinetics into 
the ground-state mode. Furthermore, the ground-state 
mode P|| = refers to the lowest energy E — and to 
the largest occupation number Ne=o- For well-developed 
BE statistics of QW excitons one has Ne=o ^ 1. As 
a result, population of the state E — generally re- 
quires an additional time in comparison with that for the 
high-energy states E > k B T with N E < 1. Thus LA- 
phonon-assisted occupation of the ground-state mode is 
the slowest, "bottleneck", relaxation process, which de- 
termines the thermalization kinetics of quasi-equilibrium 
QW excitons. A similar picture holds for thermalization 
of excitons in bulk semiconductorsoti 

With the substitution of N E = N E q given by Eq. (g) 
in Eq. ([?]) we derive 
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Equation (|ll]) describes the thermalization dynamics 
T = T(t) of QW excitons from the effective tempera- 
ture Ti — T(t = 0) to the bath temperature Tb. A finite 
lifetime r' of excitons, due to radiative and nonradiative 
recombination, can be incorporated in Eq. ( jlT| ) by the 
degeneracy temperature T (t) oc p2D(t = 0) exp(— t/r'). 
In next section we apply Eq. ([ll]) to the thermalization 
kinetics of QW excitons in the absence of their recombi- 
nation (V — > oo). This analysis refers to the case when 
the characteristic thermalization time Tth "C t'. In sec- 
tion IV we will use Eq. (^Tj) in order to develop a theory of 
resonant photoluminescence of quantum degenerate ex- 
citons in perfect QWs, when r' = r op t is determined by 
the intrinsic radiative lifetime tr of QW excitons. 



III. THERMALIZATION KINETICS OF 
QUASI-EQUILIBRIUM QW EXCITONS 

In this section the thermalization dynamics of QW ex- 
citons is analyzed within Eq. ( |Tl| ) for p2D = const. First 
we linearize Eq. (|ll|) about some bath temperature Tb 
to study the exponential relaxation of QW excitons with 
the effective temperature \T — T \ <C Tb- Linearization is 
not appropriate for the thermalization kinetics at T = 0. 
Here we show that the thermalization becomes very slow, 
with T(t) oc 1/lnf. In the last subsection, the cooling of 
hot QW excitons with T — > Xj, is treated for both 
classical and strongly degenerate limits of BE statistics. 
There are four characteristic temperature scales, T, To, 
Th, and Eo/ks, and the precise form of the relaxation 
will depend on all four. Note however that if T ^> To, 
the degeneracy temperature is not a relevant parameter, 
and the most relevant aspects of phonon bottleneck ef- 
fects are captured by the ratio Eo/kBT . In practice the 
relaxational dynamics is controlled by two parameters, 
Tq/T and E /k B T b . 



A. Relaxation kinetics between nearby 
thermodynamic states 

If the effective temperature T of QW excitons is close 
to the bath Tb, the basic thermodynamic Eq. ( pi]) reduces 
to 
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(T - T b ) , 



(12) 



dt T t h 

where the effective thermalization time Tth is given by 



[a^J e{e - 1) 
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F, 



1 



( e eE /k B T b + e -T /T b _ 



^Eo/k B T b 



l e eE /k B T b _ i\ 



(13) 



The linear approximation of Eq. by Eq. ( |l2|) is valid 
for \T — T b \/T b < k B T b /E and can be done for any 
Tb > 0. Equation ( [l2| ) corresponds to the exponential 
thermalization law T(t) = Tb + AT exp(—t / r t h) , where 
AT = Ti-T b and T t = T(t = 0) is the initial temperature 
of QW excitons. The thermalization time Tth uniquely 
describes all quasi-equilibrium relaxation processes in a 
gas of quasi-2D excitons. For example, the relaxation 
kinetics into the ground-state mode E — is given by 



N E 
N 



o(*) - + (N E= 



"E=0 



N 



f 



e=o) 



exp(— t/Tth), where 
E=0 - exp(T /Tp_- 1 and N f E=0 = exp(T /T h ) - 1. We 
will analyze Eq. (|l|) in the limit of classical and well- 
developed BE statistics of QW excitons, respectively. 
1. Classical gas of QW excitons (Tb ^> To). 
In this case Eq. (|l3|) yields 
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From Eq. ([fj) one concludes, as expected, that the char- 
acteristic thermalization time of the Maxwell-Boltzmann 
distributed QW excitons is indeed independent of their 
concentration p2D- In Fig. 2 we plot Tth = T t h{Tb) cal- 



culated by Eq. fll4|) for GaAs QW of thickness L z — 
100 A. The ratio E /k B Tb determines the high- and low- 
temperature limits of Tth- 
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FIG. 2. The thermalization time Tth = Ttfe(T&) of a classi- 
cal gas of QW excitons. GaAs QW with L z = 100 A and r sc 
= 13.9 ns [see Eq. (H)]. Inset: the high-temperature thermal- 
ization time th = th(L z ) oc of Maxwell-Boltzmann dis- 
tributed QW excitons (solid line) and the high-temperature 
thermalization time toh = tqh(L z ) oc L z of strongly statisti- 
cally degenerate QW excitons ( dot-dashed line). The latter 
dependence refers to T = {k B T^)/E Q = 1 (T > T b ). 

A. High-temperature limit (k B Tb 3> The char- 

acteristic thermalization time, which in this case we will 
designate by th, is given by 

1 C2D , 

— = , where 

th t sc 

2 



Q 



2D 



2tt 



de e 



1 



1 



(15) 



In the high-temperature limit of a classical gas of QW 
excitons the thermalization time t b is independent of 
the bath temperature Tb. The constant C2D = C2d(cl), 
which completely determines th in terms of r sc , is much 
larger than unity. For example, for GaAs QW with 
L z = 100 A the dimensionless parameter a = 0.096 and 
the constant C2D = 2530. The corresponding thermal- 
ization time of the QW excitons is th = 5.5 ps. For com- 
parison, the high-temperature thermalization time of a 



classical gas of excitons in a bulk semiconductor is given 
(3/87r)rjI 



u ( 3£ >) 



This estimate yields t^ D ^ 



1.66 ns for bulk GaAs, i.e., r 



ratio t 



(2D), (3D) _ 



(3D) 
H 



» ^ 



H 



(87r)/(3C 2Z5 ) 



th ■ The 
3.3 x 10" 3 , which 

refers to bulk GaAs and GaAs QW with L z = 100 A, 
clearly demonstrates effective cooling of QW excitons in 
the presence of a bath of bulk phonons. Due to the 
relaxation of momentum conservation in z-direction a 
ground-state QW exciton couples to the continuum states 
E > Eq rather than to the single energy state E = Eq 
as occurs in bulk materials. With decreasing bath tem- 
perature Tb the effective thermalization time T t h of QW 
excitons monotonously increases starting from its high- 
temperature limit given by th (see Fig. 2). For Tb = 5 K 
one finds from Eq. ( |l5| ) that Tth = 19-1 V s m GaAs QW 
with L z = 100 A. 

There is an uncertainty in values of the deformation 
potential D of exciton - LA-phowm .interaction, pub- 
lished in literature (see, e.g., Refs.all3'Ell), from D m i n = 
7.0 eV to D max — 18.1 eV. In our numerical evalu- 
ations we use D = 15.5 eV. According to Eq. Jl2|), 
the deformation potential contributes to T t h only through 
t sc oc 1/D 2 . Therefore the numerical calculations of T p h 
can be straightforwardly re-scaled to another value of D 
(we will use this procedure in calculations of Fig. 9apto 
reproduce qualitatively the experimental data of Ref.a). 
The value D = 15.5 eV corresponds to t sc — 13.9 ns. 

B. Low-temperature limit (ksTb < Eq). In this case 
Eq. (pi) reduces to 



1 

TL 



2tt 3 / 2 



En 



k B T b 



1/2 



^—Eo/ksTb 



(16) 



where we designate the low-temperature limit of Tth by 
Tl. According to Eq. (|l^) the thermalization time tl 
increases exponentially with decreasing bath tempera- 
ture Tb below Eo/ks- The origin of this result can 
be understood from the initial Eq. (^) of the phonon- 
assisted relaxation kinetics into the ground-state mode 
E = 0. For a classical gas of QW excitons one has 
Ne C 1. As a result, both the spontaneous (indepen- 
dent of Ne=o) and the stimulated (proportional to Ne=o) 
processes contribute to population of the ground-state 



mode. The first, Stokes, process is oc Ne>e o 0- ' 



,) 



and increases Ne=q, while the second, anti-Stokes, one 
is cx — -/V£ = o n B>£; and decreases the occupation number 
of the mode E = 0. For the high-temperature limit of 
a classical gas of QW excitons both opposite fluxes are 



intense because 



, ph 
l E> 



E = k B T b /E > 1. In the low- 



temperature limit n p E l >Eo = exp(— E/ksTb) <C 1 and 
both the spontaneous process, which is proportional to 
NE>E (l + n p E h >Eo ) ~ N E >E = (T /T)ex V (-E/k B T), 
and the stimulated process, which is proportional to 
-N E =Qn p E > Eo = -(T /T)exp(-E/k B T b ), are exponen- 
tially weak. 

The plot Tth = T t h{Tb) (see Fig. 2) shows that for exci- 
tons in GaAs QWs the high-temperature limit given by 
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Eq. (JT5J) means T b > 30 - 50 K. On the other hand, a 
strong increase of r t h occurs already at Th ~ 5 — 10 K, 
i.e., at temperatures much above the low-temperature 
range determined by X& < Eo/kb ~ 0.54 if and where 
the approximation of Tth by Eq. ( fL6| ) is valid. 

2. Quantum gas ofQW excitons (To 3> T{J. 

In this case Eq. (n3h reduces to 



716/, V T sc/ \KB-loJ Jl 



de e 2 



F, 



(ay/e(e-\)) 



e-l 

; ge-E /fe B T(, 
( e eE /k B T b _ 1^2 



(17) 



Equation ( |17| ) shows that the thermalization time of 
quantum degenerate quasi-2D excitons depends on the 
concentration p 2 o, i- e -, T th oc Pid- 

A. High-temperature limit (ksTb ^ Eq). In this limit 
Eq. @ yields for t oh = T th : 



1 

toh 



k B T b 2 \ C 2D 
T Q Eo l t sc 



C 2D = 2tt I d: 



where 
F z (a^e(e-1) 



e- 1 

From Eqs. (§J|) and (Q we find 

1 _ / ^bTJ 2 \ / C2D \ 1 



(18) 



(19) 



The constant C 2 d = C 2 d{o) is much larger than unity, 
but much less than C 2 u- For example, C 2 d = 170 
so that C 2D /C 2D = 0.07 for GaAs QW with L z = 
100 A. According to Eq. ( [l9| ) one has t h < t~h for 
T » T b > [(r E C 2D )/(fc s (7 2D )] 1 / 2 and t „ > t h 
for T b < [(ToEoC 2D )/{k B C 2D )}^ 2 . The acceleration or 
slowing down of the relaxation kinetics in comparison 
with that in a classical gas of QW excitons originates 
from BE statistics. Because for the quasi-equilibrium de- 
generate QW excitons Ne=o = exp(To/T) ^> Ne>e = 
l/[exp(E /kBT) — 1], only the stimulated processes cx 
Ne=o{Ne>e (i ~ n p E>E ) contribute to occupation kinet- 
ics of the ground-state mode E = 0. From Eqs. ([?]) 
and (|) one obtains ST/T = -SN E =o/N E =o for clas- 
sical statistics and ST/T = -(T/T )(8N e =o/Ne=o) for 
well-developed BE statistics of the ground-state mode, 
i.e., the same relative change of the occupation num- 
ber Ne=o is accompanied at T <C To by the much less 
relative change of the effective temperature than that 
at T To- Therefore, BE occupation of the ground- 
state mode with Ne=o 3> 1 slows down the thermal- 
ization kinetics of QW excitons. On the other hand, 
one has Ne>e 



ph 
l E>E 



<>E (ST/T b ) for well- 
developed BE statistics of the energy modes E > Eo 
and N e >e ~ uV e>e q = nP E>E ( E / k BTb)(ST/T b ) for a 
classical distribution of QW excitons at E > Eo- As a 



result, the BE occupation numbers Ne>e ^> 1 enhance 
the relaxation dynamics by the factor ~ k B Tb/ Eo, which 
is much larger than unity in the high-temperature limit. 
With decreasing bath temperature Tb the slowing down 
of thermalization, which results from N B =o 3> 1j starts 
to dominate over the acceleration of relaxation due to BE 
statistics of the modes E > Eo- 
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Tth{Tb) 

BE-distributed QW excitons of the densities 10 12 cm 
(solid line), 5 x 10 11 cm" 2 (long dashed line), 10 11 cm -2 
(dashed line), and 10 10 cm -2 (dotted line). The calcula- 
tions with Eq. (|l|) refer to GaAs QW with L z = 100 A and 
T ac = 13.9 ns. 



The thermalization time Tth — T t h{a) is very sensi- 
tive to the dimensionless parameter a cx L z through the 
form- factor function F z on the r.h.s. of Eq. (p^). With 
decreasing L z the spectral width of F z (ayje(e — 1)) in- 
creases indicating a stronger relaxation of the momentum 
conservation in the z-direction. From Eqs. (|l^) and (|Tg| ) 
we conclude that for the infinite square QW confinement 
potential t b oc L 2 z and toh oc L z , respectively. The de- 
pendences th — th{L z ) and tqh = tqh{L z ) are plotted 
in the inset of Fig. 2. 

B. Low-temperature limit {k B Tb < Eq). In this case 
we get from Eq. (|l7|) : 



1 

T0L 



2^ 3 / 2 / E T b 



1/2 



\k B T 2 



-E /k B T b 




(20) 



where we designate the low-temperature thermalization 
time of degenerate QW excitons by tql and the cor- 
responding thermalization time of Maxwell-Boltzmann 
distributed QW excitons tl is given by Eq. The 
slowing down of the thermalization process by the fac- 
tor TojT b oc Pid, which is much larger than unity, 
stems from well-developed BE statistics of the ground 
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state mode, i.e., from N B =o 3> 1. The distribution 
function of QW excitons at E > Eq is classical, i.e., 
Ne>e = exp(—E /k B T) < 1, and does not enhance the 
relaxation processes into the ground-state mode E = 0. 
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FIG. 4. The thermalization time rth ~ Tth{p2D) of 
BE-distributed QW excitons at the bath temperatures 3 K 
(solid line), 5 K (dot-dashed line), 10 K (long dashed line), 
20 K (dashed line), and 50 K (dotted line). The calculations 
with Eq. dij|) deal with GaAs QW with L z = 100 A and 
t sc = 13.9 ns. Note that numerical evaluations of Figs. 2-4 
refer to effective temperatures of QW excitons close to the 
bath temperature, i.e., \T — Tb\ -C T&. 

The influence of BE statistics on the thermalization of 
QW excitons is demonstrated in Fig. 3 for r t h = Tth{Tb), 
P2D = const., and in Fig. 4 for r th = T th (p 2 D), T b = 
const., respectively. The numerical evaluations of r t h 
refer to Eq. @. 



B. Thermalization of QW excitons at Tj, = 



The linearization of the basic Eq. (11) at T& > 0, 
developed in previous subsection, does not hold when 
\T — Tb\ > k B T 2 /Eq and especially at zero bath tem- 
perature. Because N E q =Q — > oo at Tb = 0, the effective 
temperature T of QW excitons is a nonanalytic point of 
the r.h.s. of Eq. Q at T = T b = 0. For T b = Eq. @ 
reduces to 



d_ 

dt' 



T : 



-(t 

E 



-To IT 



i V e- 1 



(ay/e{s 



e eE /k B T . 



-T /T _ I 



(21) 



Equation (^l]), which describes how the QW excitons 
with effective temperature T cool down towards Tb = 0, 
can be further simplified for ksT < Eq : 



dt 



-T 



T 2 



k B T 
E Q 



1/2 



-E /k B T 



(22) 



In dimensionless time and temperature units r = 
(2tt 3 / 2 E t)/{k B T T sc ) and f = k B T/E Eq. (f§) takes 
the canonical form cif/rfr = -f 5 / 2 exp(-l/f). The so- 
lution T — T(t) of Eq. ((H) is given by the transcendental 
equation: 



F{k B T/E ) = 



2tt 3 / 2 



En 



k B T 



t + Ai 



(23) 



where F(x) = [1/y/x - Ds{l/ y/x)]e 1/x , Ds(y) = 
dt e is Dawson's integral, and the integration 

j 

constant Ai is defined by the initial condition Ti — 
T(t = 0) < E a /k B , i.e., A, = F{k B TjE a ). For 
t ^> JYA:BTo)/(27r 3 / 2 £ ; o)]T sc the asymptotic solution of 
Eq. (E3) is 



k B T(t) 



En 



In 



(2^ 3 / 2 £;ot)/(fcBToT sc ) 



(24) 



i.e., f (r) = l/ln(r) for r > 1. 

The asymptotic law Eq. ( p4] ) characterizes the nonex- 
ponential and extremely slow thermalization kinetics of 
a quantum gas of QW excitons. Because the phonon 
occupation numbers n^ 1 = at Tb = 0, only Stokes scat- 
tering of the QW excitons determines the relaxation pro- 
cess. In this case the integrand on the r.h.s. of Eq. (||) 
is proportional to Ne>e (Ne=o + 1) — Ne>e Ne=o- 
At k B T < En the energy states E > Eq, which cou- 
ple with the ground-state mode E = [see Eq. (||)], are 
weakly populated by the QW excitons, while ^^=0 ^ 1- 
The Maxwcll-Boltzmann distribution function Ne>e = 
cxp(—E/k B T) gives rise to the factor ex-p(—En/k B T) <C 
1 on the r.h.s. of Eq. j22]). This term together with the 
need to accumulate a huge number of QW excitons in 
the ground-state mode pn = 0, due to BE statistics, are 
responsible for l/ln(r) thermalization law at Tb = 0. 

In Fig. 5 we plot the results of numerical evaluation of 
the relaxational dynamics at Tb = from Eq. ( pl[ ) for var- 
ious concentrations p 2 D of QW excitons and for an initial 
effective temperature Tj = 100 K. The simulations re- 
fer to the long-lived quasi-2D excitons in GaAs/AlGaAs 
CQWs. The first hot relaxation (see the inset of Fig. 5) 
is completed within At < 30 — 100 ps. At the end of the 
transient stage the effective temperature T of QW exci- 
tons is still much higher than Enjk B ~ 0.54 K so that 
Ne~e > 1 for T < Tq. As a result, the duration of hot 
thermalization decreases with increasing concentration 
Pid of QW excitons. Here one has an acceleration of the 
relaxation kinetics due to BE statistics of energy modes 
E > Eq. The further thermalization kinetics of cold QW 
excitons refers to t > 100 ps and reveals 1/ ln(r) law (see 
Fig. 5). According to Eq. (p4j), critical slowing down of 
the phonon- assisted relaxation dynamics develops with 



8 



increasing degeneracy temperature Tq oc p2D- An exten- 
sion of the numerical evaluations presented in Fig. 5 on 
ps time scale shows that even at At — 1 ps the effective 
temperature of QW excitons is still a few hundred mK for 



10 1 



T(t 



1/xs) = 89 mK for piD = 

ill „™-2 



P2D > 

10 10 cm-' 2 , T(t = 1/is) = 129 mK for paD = 10^ cm 
T(t = lfis) = 181 mK for p 2D = 5 x 10 11 cm~ 2 , and 
T(f = 1/xs) = 213 mK for p 2 D = 10 12 cmr 2 . On a 
time scale much longer than the duration of the initial 
hot thermalization the effective temperature T is nearly 
independent of T t = T(t = 0), provided that k B Ti E . 
Recently, a strong increase of the thermalization times 
in a highly degenerate quantum gas of quasi-2D excitons 
has indeed been observed in GaAs/AlGaAs CQWsc3. 
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FIG. 5. Thermalization dynamics T = Tit) of QW exci- 



tons at T), = 0, p2D 



10 1 



(solid line), 5 x 10 1 



(long dashed line), 10 cm (dashed line), and 10 cm 
(dotted line). Inset: the initial transient thermalization of 
hot QW excitons BE-distributed at t = with Ti = 100 K. 

However, our treatment has neglected possible low- 
temperature collective states due to interactions between 
QW excitons. At a critical temperature T c = aTg, where 

(on") 1/2 

a = a(a x p 2D ) < F a system of quasi-2D excitons may 
undergo a phase transition to a superfluid state. Our cal- 
culations of relaxational thermodynamics become invalid 
at T < T c , i.e., one cannot trace with Eq. ( pl| ) the tran- 
sition to the excitonic superfluid phase and how further 
a collective ground state of QW excitons arises at T — > 
(for P2d[o!x D ^] 2 < 1 the ground state can be-iaterpreted 
in term of BE condensation of QW excitonscZEa) . While 

in the dilute limit ^(p2d[&x' D ' > ] 2 ) "C 1 the parameter 
a«l, i.e., T c -C 7ot3, there is still no first-principle the- 
ory to estimate a for the densities P2d[ox 2D " > ] 2 < F With 
the mentioned above restrictions, we find that cooling of 
high-density QW excitons to very low temperatures is 
rather difficult. This conclusion may have some bearing 
on the searcbfc!3 for the collective ground state of indirect 
excitons in GaAs/AlGaAs CQWs. 
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FIG. 6. Fhe relaxational dynamics of a classical gas 
of hot QW excitons at T b = 1 K (solid line), 4.2 K 
(dot-dashed line), 10 K (long dashed line), and 25 K (dashed 
line). Fhe latter straight line refers to the exponential ki- 
netics In |(T - T b )/(Ti -Tb)\ = -t/TH which is valid in the 
high-temperature limit ksT^ 3> Eq. 



C. Thermalization dynamics of hot QW excitons 
(T > T b ) 

In experiments [1-9] with nonrcsonant excitation of 
heavy-hole excitons in GaAs QWs the initial effective 
temperature Tj 3> TJ,. While in this case the thermaliza- 
tion kinetics can be analyzed numerically from the basic 
Eq. (|ll]), we will clarify the various relaxation scenarios 
by analytic approximations of Eq. (pd|). 

1. Classical gas of QW excitons (Ti 7], 3> Tq). 

A. High-temperature limit (ksTb S> Eq). In this case 
Eq. (O) reduces to 



d 1 , 

dt T = -^tT-n) 



(25) 



where the characteristic thermalization time th is given 
by Eq. (|i"5|). Thus, in the high-temperature limit ksTb ^> 
Eq the exponential law T(t) = T b + (Ti ~ T b ) exp(-t/r H ) 
is valid even for (T* - T b )/T b > 1. 

B. Low-temperature limit (ksTb < Eq). The initial 
hot thermalization down to fc^T ~ Eq is approximated 
by 



d T _ T 
dt t h 



(26) 



This exponential cooling T(t) = Ti exp(— t/rjj ) com- 
pletes at t' ~ th ln(Ti/Tb). The hot thermalization refers 
to the high-temperature limit given by Eq. (p5|). At t > t' 
the cold QW excitons relax to the bath temperature ac- 
cording to 



^-T = --(T-Tb) 
at t l 



(27) 
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where the low-t emp erature thermalization time tl is 
defined by Eq. (16). Approximation (p7|) is valid for 
(T — T b )/T b < k B T b /Eo, yields the exponential relax- 
ation with tl, and corresponds to the linearization of 
Eq. (0) at T = T b . 

If the high-temperature limit does not hold, the ther- 
malization kinetics of Maxwell-Boltzmann distributed 
QW excitons is exponential only locally in time, i.e., 
i~th = i~ t h{t). The described above picture of the relax- 
ation in the low-temperature limit in fact deals with a 
continuous increase of the instant thermalization time 
Tth{t) from tjj to tl- As a result, for Tj — T b ^> T b the 
total relaxation process can show a nonexponential be- 
havior. In Fig. 6 we illustrate how a single-exponential 
thermalization kinetics with Tth = t~h of a classical gas 
of QW excitons initially distributed at T T b devel- 
ops with the decreasing bath temperature towards the 
nonexponential relaxation. 
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FIG. 7. Thermalization kinetics of hot QW excitons of 
the densities 10 12 cm~ 2 (solid line), 5 x 10 11 cm -2 (long 
dashed line), 10 10 cm -2 (dashed line), and 10 9 cm -2 (dot- 
ted line). The initial effective temperature of QW excitons 
Ti = 50 K, the bath temperature T b = 4.2 K. GaAs QW 
with L z — 100 A and r sc = 13.9 ns. The first transient 
thermalization occurs in the time interval I (t < 20 — 40ps). 
The density-dependent quasi-steady-state relaxation develops 
in the time-domain //, where T — Tb < Ti,. Inset: verification 
of the approximation Eq. (^9|) of nonexponential transient re- 
laxation of strongly degenerate QW excitons. 

2. Quantum gas of QW excitons (Ti Tq^$> Tb). 

A. High-temperature limit (k B Tb ^> Eq). The initial 
hot relaxation to the temperature T ~ To is given by 
Eq. (fH and completes at t' ~ t h ln(T/T ). At t > t' 
QW excitons become degenerate and BE statistics starts 
to influence the thermalization process. The relaxation 
kinetics is approximated by 



where t h oc p2D is given by Eq. (|18|). Equation ( |2§| ) 
yields the nonexponential transient thermalization from 
T ~ T to T ~ T b : 



T(t) 



l + 2{T /T b Y[{t-t')/T QH 



1/2 



(29) 



^T=-(L 
dt \T b 



(T - T b ) 

Tqh 



(28) 



The nonexponential transient lasts for a time At ~ 
tqh/2. The relaxation kinetics from t > t" ~ t' + Tqh/2 
refers to T — T b < T b and recovers the exponential law 
with tqh- 

B. Low-temperature limit (fc^Tb < Eq)- At < 
t < t' ~ r/fln(Tj/To) the first hot thermalization 
from T = T to T ~ Tb is characterized by the ex- 
ponential law of Eq. (p6|). The next nonexponential 
transient relaxation from T ~ To to T ~ Eo/ks is 
given by Eq. (K9) and corresponds to the time inter- 
val t' < t < V ~ (k B T b /E ) 2 (T 0H /2). These two 
stages are similar to those found in the high-temperature 
limit. For the thermalization from T ~ Eo/ks to 
T ~ T b one recovers Eq. $22\) and the corresponding 
1 / ln(r) kinetics of Eq. (|24|) treated in the previous sub- 
section for T b = 0. This stage is rather long and 
completes at t'" ~ (k B T T sc /27r 3 / 2 E )exip(Eo/k B T b ) = 
(ksTh/ ' EqY^tql- The final exponential relaxation refers 
to (T - T b )/T b < k B T b /E and is given by T(t > 
If") = T b + (k B T 2 /E )exp[~{t - t'")/T 0L l where the 
low-temperature thermalization time tql of the degen- 
erate QW excitons is defined by Eq. (|20|). 

In order to demonstrate how BE statistics can influence 
the thermalization kinetics of QW excitons at T ^ T b , in 
Fig. 7 we plot the numerical solution T = T(t) of Eq. Q 
for Tj = 50 K, T b = 4.2 K, and various values of the con- 
centration p2D- The development of the nonexponential 
relaxation with the increasing degeneracy temperature 
Tn oc p2D is clearly seen. The insert of Fig. 7 illustrates 
the transient thermalization given by Eq. (|29|). 



IV. RESONANT PHOTOLUMINESCENCE OF 
QUANTUM DEGENERATE QUASI-2D 
EXCITONS 

In high-quality GaAs QWs the decay of quasL2D exci- 
tons is mainly due to radiative recambinatiorjjlij. More- 
over, following the first studiesEil of the optical prop- 
erties of excitons in GaAs/AlGaAs CQWs, very recent 
experimentalj show that the radiative recombination 
channel can also be dominant for the long-lived indirect 
excitons. In this section we generalize the theorySO of 
steady-state resonant PL of a quasi-equilibrium classical 
gas of QW excitons to well-developed BE statistics and 
exploit the relaxation thermodynamics of Eq. (|ll]) in or- 
der to analyze the PL kinetics of QW excitons. 

In perfect QWs an exciton can emit a bulk photon 
only from radiative modes, which are located inside the 
photon cone given by k = k(u>) = (^fe^Lo) / (Tic) . Here 
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oj is the frequency of light and e& is the background di- 
electric constant for an exciton line. This means that 
the radiative zone of QW excitons is given by pu < ko, 
where ko = k(uj t ) corresponds to crossover of the photon 
and exciton dispersions and hut is the exciton energy. 
The intrinsic radiative rates for in-plane transverse (T- 
polarized) and in-plane longitudinal (L-polarized) dipole- 
active QW excitons are given byE3 



r T(pn) = r c 



: o P\\ 



and r L (p||) = r ( 



k 



(30) 



respectively. The intrinsic radiative lifetime tr of a QW 
exciton is defined by tr = l/r , where T = T (L Z ) — 
Tx(p\\ = 0) = Fl(p|| = 0) is the radiative rate for the 
ground-state mode p|| = 0. Z-polarized heavy-hole QW 
excitons are foxbidden. Here we use the polarization no- 
tations of Refill 

The total radiative decay rate T opt = 1 /T opt of a gas of 
quasi-equilibrium BE-distributed QW excitons is 



opt 



1 

P2D 



Ft(p\\) +Tl(p\\) 



p\\ d n 

2tt 



(31) 



where the occupation number N p ^ = is given by 

Eq. (^). Equation ( fil] ) takes into account the equal prob- 
abilities of the T- and L-polarizations and two-fold spin 
degeneracy, a + and er~, of dipole-active QW excitons. 
Using Eqs. (|) and ©) we find from Eq. ©: 




(32) 



where A = A(T, T ) = e Ek o/ k ^ T /(l - e - T °/ T ) and E ko = 
h 2 kl/2M x (for GaAs QWs with e b = 12.9 and huj t = 
1.55 eV one has E ko ~ 101 peV and E k Jk B ^ 1.17 K). 
The contribution of two-fold dipole-inactive (triplet) QW 
excitons to the total density p2D is included in Eq. ( p^ ) 
through the degeneracy temperature To oc p2D- 

For Maxwell-Boltzmann distributed QW excitons at 
the effective temperature T much larger than To and 
E k Jk B we find from Eq. © that J T = 3J L = 
{% 2 kl)/{M x T) and 



pel 
1 opt 



2 1.2 



h z k 



3M x k B T 



To 



(33) 



In the limit T <C T of a strongly degenerate gas of QW 
excitons one approximates Jt = 1 + (T/T ) ln(4E ko / kbT) 



and J L = 1 + (T/T ) ln(4E fco /fc 6 T)-2 . In this case T opt 
of Eq. (|3^ ) reduces to 



r 9 - 

L opt 



T T (^J4_ 

T To I M x k B T 



4 ■ (34) 



The expansion of Jl(T, To), Jt(T, To), and F OJ , t in terms 
of the ratio T/To <C 1 is valid for the effective tempera- 
ture of QW excitons T < E ko exp(-T /T). The second 
and third terms in the square brackets on the r.h.s. of 
Eq. (||) are small corrections, i.e., r* pt (T -> 0) -> T /2. 
The nonzero limit of T opt at T — > 0, which is com- 
pletely determined by the intrinsic radiative rate To of 
the ground-state mode p|| = 0, is due to the effective ac- 
cumulation of low-energy QW excitons with N B 3> 1 in 
the radiative zone p|| < fco- 

In the limit T 3> Tp , E kn of a classical behavior one 
derives from Eqs. ©-©: 



3M x k B T 



2 1,2 



h k. 



]TR+ To 



3M x k B Tp 
2h 2 kl 



(35) 



l/r pt we keep not 



where in the expansion of T^ pt 

only the leading term 1/Tf pt oc T [see Eq. (33)], but 



also the temperature-independent correction. The next 
terms of the expansion are proportional to i/T""-- 1 — > 
and can indeed be neglected. The first term on the r.h.s. 
of Eq. ©, i.e., l/r£L= (3M x k B Tr R )/(h 2 k^), is a well- 
known result of Ref.liS. The temperature-independent 
correction, which is given by the second term on the 
r.h.s. of Eq. (|35|), ori gina tes from the BE distribution 
function used in Eq. (31) and consists of the density- 
independent and density-dependent contributions of the 
opposite signs. The density-dependent contribution can 
be much larger than tr if k B T 3> k B To 3> E k[) . In the 
leading term 1/r^L oc T one has that the optical de- 
cay time Tgp t 3> Tr, because only a small fraction of 
QW excitons oc cup y the radiative modes p\\ < ko at 
T»T ,£ fco /fcjM 

In Fig. 8 we plot the effective decay time of quasi- 
equilibrium QW excitons r opt = T opt (T), calculated nu- 
merically with Eq. (32) for various values of the con- 
centration p2D- Following Ref.Eil, for a GaAs QW with 
L z = 100 A we use t r = 25 ps. For T > T ,E k Jk B 
the linear behavior given by Eq. ([TJ]) is clearly seen (dot- 
ted, dashed, and long dashed lines in Fig. 8). The ref- 
erence bold solid line indicates 1/T^ pt oc T law. Note 
that the linear asymptotics of T opt = T opt (T, pio) from 
the classical range T 3> To , E ko j k B down to low tem- 
peratures T — > reveal nonzero density-dependent val- 
ues. Thus one concludes that even at high temperatures 
T S> To the influence of nonclassical statistics can be 
found in the PL of QW excitons. At T — * the limit 
rf pt -> r fl [9/10 - (3fc B To)/(4T4 )] given by Eq. © 
breaks and all curves approach the double intrinsic ra- 
diative lifetime 2tr. The density-dependent deviation 
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from the linear classical law of Eq. ([ 
increasing p2D- For example, at P2D 



develops with 
= 5 x 10 11 cm~ 2 



the quantum asymptotics r| t (T — > 0) = 2tr can be 
traced already for the effective temperature T 
(see Fig. 8). 
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FIG. 8. The optical decay time T op t versus the effective 
temperature T of QW excitons: p2D = 10 8 cm -2 (dotted 
line), 5 x 10 1(1 cm -2 (dashed line), 10 11 cm -2 (dot-dashed 
line), and 5 x 10 11 cm -2 (solid line). The intrinsic radiative 
lifetime of QW excitons is tr — 25 ps. The bold straight line 
is given by 1/Tf pt = (3M x k B T)/(h 2 k 2 ). 

According to Eq. (p0|), the decay rate of T-polarized 
QW excitons diverges at pu — ► fcoE£l. In Appendix B we 
show how the homogeneous linewidth of low-energy QW 
excitons KT hom = %/t x - x + h/r ph + H/t r ~ H/t x - x re- 
moves this divergence [see Eq. (B6)]. In the presence 
of homogeneous broadening, due to dominant exciton- 
exciton scattering, Tt(p\\) of Eq. (30) changes only in a 
very narrow band of QW states given by < fco — pu < 
7^0, where 7 = 1/ (r x ^ x ujt) <C 1. In particular, Tt(jp\\ — > 
ko) = To/(27 1 / 2 ) rather than diverges. Furthermore, 
homogeneous broadening of QW excitons from the ra- 
diative zone only slightly renormalizes the PL efficiency: 
r c J pt changes on T c J pt [l - (3/4)7 1 / 2 ] (for details see Ap- 



pendix B). The correction is even less for T q opt of Eq. (34) 



Therefore we conclude that the corrections to the optical 
decay rate T opt due to the homogeneous linewidth T; lom 
can indeed be neglected. Note, however, that inhomoge- 
neous broadening strongly influences the PL process and 
changes our results obtained for a perfect QW. 

The optical decay of QW excitons from the radiative 
zone pu < fco does not violate the assumptions of relax- 
ational thermodynamics provided that t x _ x <C Tr. In 
this case the thermal quasi-equilibrium distribution of 
the finite lifetime QW excitons at E < Ek holds due 
to exciton-exciton scattering. Because both the minimal 
thermalization time Tth given by tjj of Eq. (|l^) and the 
intrinsic radiative time 77? are on the same time scale 
of 5 — 50 ps, the condition t x _ x -C th [see Eq. (0)] 



of the thermodynamic picture also nearly guarantees 

T x - X "C TR. 

The radiative decay leads to the continuous decrease of 
the density of QW excitons (or the degeneracy tempera- 
ture To oc P2d) and gives rise to the PL signal according 
to 



° _ , , _ P2D 

tt:P2D — -1 opt{J-0, 1 P2D = -, ^7 , 

Ot T op t{P2D,T) 

T % dp 2 D _ % P2D 

J- pl = —n^>t — 7^ — = hut 



at 



T op t{p2D,T) 



(36a) 
(36b) 



= IpL(t) is the phot ol umin escence intensity. 
(|l|), (ph, and J36a|)-(|36b^ describe the PL 



where Ipl 
Thus Eqs. 

kinetics of QW excitons cooling from the initial effective 
temperature Tj to the bath T},. For a high-temperature 
classical gas of QW excitons the optical decay time rf pt 
of Eq. (pq) is approximated by the density- inde pendent 



1 /Tf pt of Eq. ( |33|) and the solution of Eq. ( B6a ) is sim- 



ply p2D{t) = p^exp(- $*Tf pt {t)dt), where p^ is the 
initial concentration of QW excitons and T^ pt is time- 
dependent through the effective temperature of QWex- 
citons T = T(t). In the general case, however, Eq. (36a) 
gives a nonexponential decay of the density. The char- 
acteristic times of the both fundamental processes which 
contribute to the PL kinetics, LA-phonon- assisted ther- 
malization and optical decay, generally depend upon the 
effective temperature and density of QW excitons, i.e., 
nh = T t h{T,T ) and T opt = T opt (r, T ). As a result, the 
PL kinetics given by Eqs. @, (||), and (]36a|)- (|36b|) is 
also T- and p2D-dependent. The thermalization process 
and radiative decay work in the opposite directions with 
respect to the quantum-statistical effects. Cooling of QW 
excitons is accompanied by an increasing number of low- 
energy particles towards well-developed BE statistics. In 
the meantime T op t decreases and the optical decay of P2D 
speeds up resulting in decreasing To- Thus the radiative 
processes interfere with the development of BE statistics. 

For a classical gas of QW excitons in the high- 
temperature limit fc gTj, ^> E q the PL dynamics of 
Eqs. (|Ti"|), ([32]), and (36a)-(36t) can be analyzed analyt- 
ically. In this case, T{t) = T b + (T { - T b ) exp(-t/T H ), 
where the temperature-independent th is given by 
Eq. (|l5|). At the first transient thermalization of hot QW 
excitons from the initial distribution at effective Ti 3> T b 
one has that rf pt {Ti) 3> t#, i.e., the optical decay is very 
slow and practically does not change the concentration 
of QW excitons. The transient stage lasts a few th and 
for this time domain Eq. (|36h|) reduces to 



(ft) 



i PL {t < t tr ) = ho; t ( o 5fL 1 pfy'™ 



hujt 



v 3M,fc B T v 
(o) 

P 2D gt/T 
T opt(Ti) 



(37) 



i.e., IpL(t < ttr) oc exp(t/T#). Therefore at the hot 
relaxation stage t < t tr the PL intensity exponentially 
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increases with the rate 1/th, due to the population of 
the radiative zone. The transient stage completes when 
T(ttr) — Tb < Tb. The PL intensity gets its maximum at 
t = to > ttr, when the population efficiency of the modes 
P\\ < fco is compensated by the increasing optical decay: 



t = t h m — I = th In 



TH 



'opt 



(Ti) 



(38) 



The corresponding effective temperature of QW excitons 
is T(to) - T b [l + T H M l pt (Ti)], i.e., T(t ) - T b « T b . 
At i = t s > to the PL dynamics becomes steady-state 
and is determined by the optical decay, because the QW 
excitons are already thermalized at T = Tb. For this time 
domain the PL kinetics is given by 



I PL (t > t s ) = hw t \ 



where = /9 2 d(*«) 



h 2 k 2 



3M x k B T b 



(s) 
P2D 



opt 



(T b ) 



(39) 



(o) 



Thus at the steady-state 
regime I PL (t > t s ) cx exp[-i/r^ t (T 6 )], i.e., the PL in- 
tensity decreases exponentially with the optical decay 
rate T^ pt (T b ). Note that because usually Tj 3> T and 
To(t — > 00) — > 0, the very first transient stage and the 
very last steady-state decay of excitonic PL always follow 
Eqs. ( |3"7| ) and (|3^), respectively. 

In Fig. 9a we plot numerical simulati ons of the PL 
dynamics done with Eqs. (pi]), (|32|), and (36a)- (36b) for 
excitons of various initial densities in GaAs QW with 
L z = 100 A. In order to model the experimental data 
on the density-dependent PL kinetics plotted in Fig. 2 of 
Ref. [4a] we use r sc = 85 ns and tr = 40 ps. The change 

of the PL dynamics with increasing p^D shown in Fig. 9a 
reproduces qualitatively the corresponding experimental 
observations of Ref. [4a]. The decrease of the rise time 
to of the PL signal with the increasing initial density 
of QW excitons is mainly due to increase of the optical 
decay rate at high densities (see the inset of Fig. 9a). 
This effect originates from BE statistics of QW excitons. 

In Fig. 9b we show numerical modeling of PL kinet- 
ics in GaAs/AlGaAs CQWs with the intrinsic radiative 
lifetime of long-lived excitons tr = 30 ns. Here BE 
statistics strongly influences the PL temporal behavior 
in the steady-state regime at t > t s , but still before a 
time when indirect excitons become Maxwell-Boltzmann 
distributed due to continuous decrease of P2D- At this 
time domain the thermalization kinetics of high-density 
CQW excitons undergoes critical slowing down, the ef- 
fective temperature of CQW excitons is nearly stabilized 
at To > T > T b , and the occupation numbers of the 
radiative modes N Pll <k a S> 1- As a result, with the in- 
creasing initial concentration the decay rate of the 
PL signal approaches r^ pt — 2tr, i.e., one has an acceler- 
ation of the PL optical decay owing to well-developed BE 



statistics [see Fig. 9b, where the long dashed and solid 
reference lines indicate ln(/pi/Jo) = (to — t)/2T opt (Tb) 
and \n(I PL /I ) = (t - t)/2r R , respectively]. 
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FIG. 9. Phototuminescence dyn amic s of Q W excitons 
modeled with Eqs. @, @, and (p6a|)-(p6b|). (a) GaAs 
QW with L z = 100 A, D = 6.3 eV, r sc = 85 ns, 
T b = 5 K , P2D = 10 9 cm -2 (dashed line), 
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Tr = 40 ps, l b = t) ft, p2D = W cm 
10 11 cm -2 (dot-dashed line), and 5 x 10 11 on' 2 (solid line). 
Inset: change of the optical decay time T op t = T op t(t). (b) 
GaAs/AlGaAs CQWs with r sc = 13.9 ns, r R = 30 ns, 
Tj, = 3 K, p2D = 10 9 cm~ 2 (bold dashed line), and 
5 x 10 11 cm~ (bold solid line). The long dashed and 
solid straight lines refer to the exponential kinetics given by 
I PL {t) oc exp[-t/T OI ,t(Tb)] and I PL (t) oc exp[— t/(2T R )], re- 
spectively. The PL intensities IpL(t) of (a) and (b) are nor- 
malized by Jo = 248.31 kW/cm 2 . 
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V. DISCUSSION 

There are three neglected factors, the interface polari- 
ton effect, the low-energy biexciton states, and residual 
interface disorder, which can affect the developed model. 

In perfect QWs, the radiative modes pj| < k(uf), which 
are responsible for the optical decay of QW excitons 
in bulk photons, are accompanied by confined modes 
P|| > k(u) located outside the photon cone. These modes 
give rise to QW interface polaritons. While the QW 
polaritons cannot be seen in standard optical experi- 
ments, which deal with bulk photons, they can modify 
the quadratic in-plane dispersion of QW excitons and 
contribute to the total optics in a "hidden" way. There- 
fore in the general case both conjugated phenomena, the 
radiative decay into bulk photons and tliCpQW polari- 
ton effect, should be treated simultaneously^. The dis- 
persion law of T-polarized QW polaritons is c 2 p 2 /cb — 

iT x uj - 



GaAs QW with L z = 100 A the interface polariton effect 
plays no role at concentrations p2D > 3 x 10 



10 



cm 



1 



[Rc^p] - e b cj 2 /c 2 )/(u; 2 + T^ t p\/M x 

to 2 ) , where R c = (ire 2 f xy ) / '(e&mo), f xy is the oscillator 
strength of exciton-photon coupling per unit area of a 
QW (the subscript xy refers to the in-plane polarization 
of the light), and ^-iM the rate of incoherent scatter- 
ing of QW excitonsoEix Actually the dispersion equa- 
tion can be applied to the both radiative and confined 
modescj. Within the photon cone the solution lo = Qj(p\\) 
is u = ujt — i^T(p\\)/2, where Tt{p\\) is given by Eq. ( J30| ) 
with To = ( v / eb/c)i? c . The dispersion u> — u>t{p\\) of 
T-polarized QW polaritons is located outside the pho- 
ton cone, approaches the quadratic dispersion %p 2 /2M x 
at p|i 3> fcoj an d continuously transforms to the photon 
dispersion cpu/ ' ■ s fth at < kg. Thus during thermaliza- 
tion the QW excitons can accumulate at the "bottleneck" 
band of the T-polariton dispersion rather than in the ra- 
diative zone given by u> — u t + Tip 2 /2M X for p\\ < fco- 

The polariton effect will not influence the thermaliza- 
tion and PL kinetics provided that E ka > (l/2)htt® w , 
where the effective QW polariton parameter Q® w is 
given by (tt® w ) 2 = {^,/c)R c bj t . In this case the QW 
polariton bottleneck band is very narrow, located at en- 
ergies above Ucot, and the mode p|| = preserves the 
status of a ground-state mode. For a CjaAs QW with 
L z = 100 A one has f xy ~ 5 x 10~ 4 A" 2 ^. This value 
of the oscillator strength yields h 2 R c ~ 1.2 x 10~ 2 eV 2 A, 
t r = l/r ~ 25 ps, and (l/2)hQQ w ~ 3 meV, i.e., 
(l/2)hnQ w > E ko ~ 0.1 meV. In perfect CQWs the 
oscillator strength f xy can be made two or three orders 



Alternatively, in the presence of the well-developed po- 
lariton picture with E ka < (l/2)hQQ w , the relaxation 
kinetics should be reformulated in terms of "QW polari- 
ton ± bulk LA-phonon <-> QW polariton" . Furthermore, 
in this case relaxational thermodynamics becomes invalid 
due to the absence of steady-state quasi-equilibrium in a 
gas of QW excitons (QW polaritons) . 

Next we consider the formation of quasi-2D biexcitons. 
The potential for exciton-exciton interaction is repulsive 
for indirect excitons in CQWs. However, in a single QW 
two excitons in some spin configurations, e.g., a + - and 
cr~ - spin-polarized, attract each other due to the ex- 
change Coulomb interaction. This interaction may lead 
to biexciton states which lie below the excitonic energies 
E > 0. While the biexciton states are absent in CQWs, 
they are presented in single GaAs QWs at low tempera- 
tures and moderate optical excitations. Typical values of 
the biexciton binding energy Em are about 1 — 2^neV 
for L z = 200 - 100 A. The law of mass actionlia ap- 
plied to QW excitons and biexcitons, which are Maxwell- 
Boltzmann distributed at the same effective temperature 
T, shows that the biexciton states are ionized in unbound 
excitons at k B T > E^ D \ i.e., at T > 10 K. At low 
temperatures, with increasing dimcnsionless parameter 
[ttm D 'f P2D the biexciton binding energy decreases due 
to screening by quasi-equilibrium QW excitons, which 
remain well-defined quasi-particles for [a x 2D ^] 2 p2D < 1. 
Here flm D ' (> a x 2D ^) is the radus of a QW biexciton. 
The screening of the biexciton states is mainly due to 
the repulsive interaction between QW excitons of the 
identical spin structure^ adiich is much stronger than 
the attractive potentiaO'Ej. The biexciton state is af- 

Em°\ where 



fected by excitons if [aE 1 *'] 2 P2dE x zu> > 
E X 2D ^ is the two-dimensional excitonic Rydberg. This 
estimate refers to low effective temperatures T — > 0. The 
squared biexciton radius is approximated by [om D 'f = 
2^ x /M x ){E x 2D) /E [2D) )[a x 2D) ] 2 ~ 2A[a x 2D) } 2 , where p. x 
is the reduced exciton mass. We estimate that the biex- 
citon state starts to weaken at 



? (2D) 



,(213) 



} 2 p 2D > 0.05-0.1, 
Certainly in the high-density 
3 x 10 12 cm" 2 , where the Mott parame- 



10 



li 



smaller than the used above value. This means that in- at concentrations p 2 D ~ 10 



i.e., at p 2 D — 
limit p2D — 2 

ter P2d[o^ 2D ^] 2 approaches unity, the thermalization and 
photoluminescence of strongly degenerate QW excitons 
occurs in the absence of the biexciton states and can be 
analyzed with our approach. During the optical decay, 

i~ 2 , the QW biexcitons 



deed E ko > (l/2)hQ9 W and one can apply Eqs. ((TTJ) , 
(|32l), and (36a)- (36b) to the relaxational and PL dynam- 
ics of indirect CQW excitons. 

The exciton-exciton interaction relaxes the polariton 
effect in perfect QWs and removes it, if T x — 1/t x _ x > 
Q® w . With the estimate of 1/t x _ x given by Eq. (A4) of 
Appendix we conclude that for direct excitons in a single 



can contribute to the relaxation and PL processes. For 
these densities of QW excitons the relaxation thermody- 
namics given by Eq. ( |ll"|) can be generalized to include 
the conversion of QW excitons to QW biexcitons and the 
LA-phonon-assisted relaxation within the both excitonic 
and biexcitonic bands. In this case the quasi-equilibrium 
concentrations of QW excitons aad biexcitons are con- 
nected by the law of mass actior£3. 
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This result can be motivated as follows. Due to the 
anti-Stokes component, only a small phase space volume 
given by p\\ < pj lax = (2kBT b )/(hv s ) contributes to the 



(2) 



The developed thermalization and photoluminescence 
kinetics refer to a perfect QW. With weak disorder, re- 
laxation can be "disorder-assisted" , whereby momentum 
and energy conservation are satisfied by LA-phonon emis- 
sion concomitant with scattering from the disorder po- 
tential. This will be relevant for relieving the phonon 
bottleneck if there are significant Fourier components of 
the disorder with wavevectors pd > 2M x v s /fo = 1.9 x 

10 5 cm -1 . For the opposite limit of strong interface and ing the next order of smallness with respect to tql oc 



two-phonon relaxation. The dependence 1/tq L ' 



rp4 

1 b 

stems from 7V P|| ^ h (^+ a + 1) oc T b 3 and fromp™^ oc T b 
on the r.h.s. of Eq. (|40[) . The two-phonon-assisted 
thermalization time Tq^ is proportional to t 2 c indicat- 



alloy disorder, the phonon-assisted relaxation kinetics of 
localized QW-excitons is modeled in Rcf.Ej and more re- 
cently in Ref.E!3 where the PL process is also studied. The 
optical decay of localized QW excitons is investigated in 
Ref.E3. Note that the number of localized QW excitons 
saturates due to Pauli blocking, if the concentration of 
the residual defect centers No is much less than p2D (m 



perfect GaAs QWs N D < 10 



10 



In this case the 



localization processes do not affect the BE kinetics of free 
QW excitons. 

At low bath temperatures T b < Eq the large thermal- 
ization times tl oc exp(— _E /fc_BT&) and Tql = (To/T^tl 
given by Eqs. ( |l6| ) and (|2^) for statistically nondegener- 
ate and strongly degenerate QW excitons, respectively, 
can be changed by two-LA-phonon-assisted relaxation 
processes "QW exciton ± bulk LA-phonon ± bulk LA- 
phonon «-> QW exciton" . While these processes are of 
the next order of smallness with respect to the consid- 
ered QW-exciton - one-phonon scattering, they can be 
dominant in the presence of drastic slowing down of one- 
phonon thermalization. The energy-momentum conser- 
vation law in two-phonon-assisted relaxation allows for a 
QW exciton with energy E < Eq to be scattered into the 
ground-state mode p|| = only by at first absorbing of a 
LA-phonon and then re-emitting of another LA-phonon. 
The rate of the anti-Stokes - Stokes two-phonon scatter- 
ing into the ground-state mode is given by 



' ill 



q,q' 



JV p „<(nJ + l) 



(£ P ii 
x 8(E P 



TlLUq 1 - 



£„„+ q ) 2 



Tiuj ph 



pii+q)^P|| 



q'+q ) 



(40) 



where E P{{ = K p%/2M x and huj^ — fiv s q are 
the energies of a QW exciton and a bulk LA- 
phonon, respectively. The matrix element of QW ex- 
citon - bulk LA-phonon interaction is M x - p h(<\) = 
\{TiD 2 q)/ (2v s pV)] 1 / 2 F z (q z L z /2), where V is the volume 
which refers to bulk LA-phonons. For a quantum gas 
of excitons Eq. (|io| ) yields the following estimate of the 
two-phonon-assisted relaxation rate: 



1 _ k B T b 



(41) 



Equations (|c| and Q show that T$(T b ) indeed 
increases with decreasing T b slowly than Toz,(7fc), he., 
the two-phonon scattering processes become dominant 
at < T b < T b c , where T§ is given by the equation 

T 0L {T b ) = T$(T b ), For example, for T b /T = 0.2 the 
crossover occurs at T b c ~ 0.06i?o//cs = 31.3 mK. Be- 
cause many-phonon-assisted relaxation of QW excitons 
at the effective temperature T < E$/kB into the ground- 
state mode occurs only through intermediate anti-Stokes 
scattering, the relaxation dynamics at T b = described 
by Eq. ( |2l| ) is universal. 

Thus our model refers to indirect excitons in perfect 
GaAs/AlGaAS CQWs and to direct excitons in single 
QWs. In the latter case, however, the concentration of 
QW excitons should be p^D > 10 11 cm~ 2 in order to sup- 
press the interface polarit on a n d bie xciton effects. The 
equations jll]), (|32]), and (36a)- (36b) provide us with an 
universal description of the thermalization and photolu- 
minescence kinetics of QW excitons within the thermo- 
dynamic picture. Note, that the concrete characteris- 
tics of QW exciton - bulk LA-phonon interaction enter 
the basic thermodynamic Eq. (|ll| ) only through the el- 
ementary scattering time r sc and the form-factor func- 
tion F z (x)- Our approach can also be applied to QW 
magneto-excitons. In this case the quantum-statistical 
effects will be considerably enhanced, because the spin 
degeneracy factor g of Eq. (Q) reduces from g — 4 to 
9=1. 



VI. CONCLUSIONS 

In this paper we have studied the influence of Bose- 
Einstein statistics on the thermalization, radiative decay, 
and photoluminescence kinetics of a gas of QW excitons. 
The numerical calculations presented in the work deal 
with perfect GaAs/AlGaAs quantum wells. The follow- 
ing conclusions summarize our results. 

(i) We have developed relaxational thermodynamics 
for QW excitons coupled to bulk thermal LA-phonons. 
The thermodynamic picture implies that the concentra- 
tion p2D of QW excitons is larger than the critical den- 
sity P2 D (for GaAs QWs with L z — 100 A we estimate 
P2 D ~ 1.2 x 10 9 cm -1 ). The relaxational thermodynam- 
ics is given by Eq. ( p"l| ) for the effective temperature T(t) 
of QW excitons. This equation describes the thermaliza- 
tion kinetics of a quasi-equilibrium gas of QW excitons 
from the initial temperature Ti to the bath T b . For a 
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quantum degenerate gas of quasi-2D excitons, when T 
is less than the degeneracy temperature To, BE statis- 
tics yields nonexponential and density-dependent ther- 
malization. In particular, for low bath temperatures 
Tb < Eo/kB ~ 1 K the relaxation processes in a degener- 
ate gas of QW excitons slow down dramatically and the 
thermalization law is given by T(t) oc 1/lni. 

(ii) The effective radiative lifetime r opt = r opt {T, To) of 
quasi-equilibrium BE-distributed QW excitons is given 
by Eq. (|32|). The temperature- independent corrections, 
given by Eq. (|35|), to the well-known classical limit r opt oc 
T show that BE statistics can be traced in the optical 
decay of a gas of QW excitons even at high effective tem- 
peratures T. Nonclassical statistics strongly influences 
the optical decay at T < Tq. In particular, at T — ► 
the lifetime r opt approaches 2tr rather than zero as in 
the classical limit T opt oc T. Here the intrinsic radiative 
lifetime tr is determined by the oscillator strength f xy 
of exciton - photon coupling in a QW. 

(iii) The PL kinetics of QW excito ns i s mod eled 
by three coupled Eqs. @, @, and (p6a])-(p6b|) for 
the effective temperature T(i), effective radiative time 
T opt {T,p 2 D) = l/r opt (T, P2d), and concentration p2D(t), 
respectively. This approach describes within the thermo- 
dynamic picture both the rise of PL, due to thermaliza- 
tion, and the decay of PL, due to the radiative recombi- 
nation of thermalized QW excitons. For a classical gas 
of QW excitons the PL kinetics depends on the effective 
temperature T, but is independent of P2D- With increas- 
ing density p 2 _D of QW excitons the quantum-statistical 
effects builds up and the PL kinetics becomes density- 
dependent. 



where Uq/S is the potential of exciton-exciton Coulombic 
interaction in a QW, S is the area. For CQWs this po- 
tential is repulsive and independent of the spin structure 
of interacting excitons owing to a rather small contribu- 
tion of the exchange interaction. For a single QW the 
potential depends upon the spin state of excitons and 
can be repulsive (e.g., for the QW excitons with an iden- 
tical spin structure) as well as attractive (e.g., for <r + - 
and a~- spin polarized excitons)cjL However, the repul- 
sive potential strongly dominateso and determines the 
equilibration dynamics of QW excitons at p2D > P2D • 

In order to evaluate a characteristic equilibration time 
t x - x we analyze within Eq. (Al) the dynamics of a small 
fluctuation <5iV (0) = 5N (t = 0) < of the ground- 
state mode population, provided that the occupation 
numbers of all other modes are given by the BE distribu- 
tion of Eq. (||). The linearization of the r.h.s. ofEq. (Al) 
with respect to 5Nq yields an exponential law for the 
decay of the fluctuation, i.e., d(5N )/dt = —8N /t x ^ x 
and 5No(t) = SN^ exp(— t/r x ^ x ). The decay time of 
the fluctuation, which we attribute to the characteristic 
equilibration time, is given by 

1 ~ ( U °2 IX ) 2 (^r) e ~ T ° /T (l-e -2WT ) 

/•oo p2tt p 2u 

x du d<p- 
Jo Jo 



e u(l-coB^) _|_ e -T /T _ l 
1 



3 tt(l+cos0) _|_ e -To/T _ l 



+ e- T o/ T - 1 



(A2) 
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APPENDIX A: EQUILIBRATION DYNAMICS OF 
QW EXCITONS 

The equilibration kinetics of QW excitons is mainly 
due to particle-particle scattering and refers to the 
Boltzmann-Uhlenbeck equation 



d N kl> =^ £ [iV k|| iV P|| (l + iV P|| _ q|| 

Pii.qn 



at 



s 



x(l+iV k +q „)- (JV k| ,+l)(N P|l +1)JV P „_ q|| JV kl+q 



x d(E p ..-q, +£ , k|| +q,i -E k 



E„ 



(Al) 



In order to derive Eq. (A2) from Eq. (Al) we put kj| = 
and introduce new dummy variable qy' = qii — Py/2. 
The dimensionless integration variable u on the r.h.s. of 
Eq. (A2) is given by u = (ftpy) 2 / \4k B M x T), <$> is the 
angle between qy' and pn. 

While the repulsive potential [/ for QW excitons 
in an identical spin state can be explicitly written in 
terms of an integral convolution of the different pair 
Coulomb potentials between the constituent electrons 
andr-holes with the excitonic wavefunctions (see, e.g., 
Rcf.H), here we give a scaling estimate of Uo. The po- 
tential Uo is determined by Uq = J U x — x (ru)dTu, where 
U x - X (r\\) is the real-space potential of exciton-exciton 
interaction and rii is the distance between two interact- 
ing QW excitons. The characteristic energy and length 
scales for U x - X (ru) are given by the two-dimensional Ry- 

dberg E x ' = 2p x e 4 /h 2 and the corresponding excitonic 

Bohr radius a x 2D ^ = Ti 2 / (2p x e 2 ), where p x is the reduced 
mass of a QW exciton. For simple model potentials 

?7a!_a,(r||) = CE X ; exp(-r||/ai 

C^( £ u 'Vr||)«p(-r||/a£ ,D) ) I 
C ~ 1, one estimates 



') and U x - X {r\\) = 
where the constant 



1G 



U a = 2irC E x 2D \a x 2D ^] 2 = vC 



(A3) 



A remarkable feature of Eq. (A3) is that while with de- 
creasing confinement in the ^-direction (the QW growth 

direction) E X 2D ^ and a x 2D ^ start to approach the cor- 
responding bulk values, their combination E x 2D \a x 2D ^] 2 
does not change. Such behavior of Uq is due to the quasi- 
two-dimensionality of QW excitons. The constant C, 
which is of the order of unity, depends on the design of 
a QW, the shape of the model potential U x - X {r\\), etc.. 
In further analysis we put C = 1. 

For a classical gas of QW excitons, when T ^ To, one 
obtains from Eqs. (A2)-(A3): 



— = Tn ("I ^ 



(A4) 



where T is given by Eq. (Q). According to Eq. (A5), the 
scattering rate due to particle-particle interaction is pro- 
portional to the concentration pio of QW excitons. Fur- 
thermore, as a signature of the quasi- two-dimensionality 
of excitons in QWs, the characteristic equilibration time 
t x - x is independent of the temperature T (for T To) 

and of the scattering length ~ a x 2D \ 

For a quantum gas of QW excitons (T < To) we find 
from Eqs. (A2)-(A3): 

The considerable increase of the scattering rate given by 
Eq. (A5) in comparison with that of Eq. (A4) is due to 
the low-energy QW excitons with ^> 1 dX E < ksT . 



APPENDIX B: HOMOGENEOUS BROADENING 
IN THE RADIATIVE DECAY OF QW EXCITONS 

The joint density of states J(p\\) for the resonant op- 
tical decay of a QW exciton with p« < ko into a bulk 
photon is given by 



J (P\\ 




dk i 



- horn 



T 2 

horn 



(Bl) 



where Tuj) x (p\\) = + h p?,/2M x and u>j(p^,k±) = 
i c l \ftb)\Jp 2 + ^1 are the dispersions of QW excitons 
and bulk photons, respectively, k± is the z-component 
of the wavevector of bulk photons, which resonantly in- 
teract with a QW exciton with in-plane wavevector p\\ , 
and Thom is the homogeneous linewidth of QW exci- 
tons. If du! 7 (pH,kx)/dkx 7^ at k± determined by 



the equation u> x (pu) — u> 7 (p||, fcj_), the integrand on the 
r.h.s. of Eq. (Bl) preserves its Lorentzian shape in 
wavevector coordinates and the joint density of states 
J(p\\) is independent of Thom- This is a regular case, 
which docs not hold for py — > ko (pii < ko). In this 
case the solution of u) x (ko) = ujj(ko,k±) is k± = and 
duj 7 (j)\\ = ko , k± ) /dk± = 0, indicating ID van Hove sin- 
gularity in the joint density of states. As a result, J(p\\) 
becomes r^ om -dependent in a close vicinity of p» = fco- 
This is the only way that the homogeneous linewidth in- 
fluences the optical decay of quasi-equilibrium excitons 
in ideal QWs. 

With substitution k± = put&ncf), where <j) 3 
[-7r/2,7r/2], Eq. (Bl) reduces to 



J (P\0 = Z2 



1 




J 


7? 




Jo 



r ir/2 



du 



7 



(B2) 

(l/4)u 4 - (S 2 /2)u 2 + (6 4 /4 + 5 2 j 2 ) 

where 7 = T hom /uj t , S = cos0 o = P\\/k = (cp||)/(e fc w t ), 
8 = sin 4>o = yl — S 2 , and the integration variable u — 
^ — 0o + sin O . Straightforward integration of Eq. (B2) 
yields the joint density of states: 

V2 p\\ 7 



J(P\\) 



Ux(P\\) 



45 2 7 2 + £ 4 



1/2 



(4(5 2 7 2 



d 4 



1/2 



6 2 



1/2 



(B3) 



For 6 3> v^^T) which is identical to the condition k*. 
pn 3> 7&0, Eq. (B3) reduces to the standard formulaES: 



J(p\\ < k 




In opposite limit 5 <C \/2~J, i.e., for pu 
yields 



(B4) 



k , Eq. (B3) 



1 



J(P\\ = k o) = 7T- 



k 



(B5) 



27T [Thorn^t) 1 / 2 

Equation (B5) shows how the homogeneous linewidth 
Thom removes the van Hove singularity at p\\ = ko- For 



5 < V27 the joint density of states is r?,, om -dependent. 

With Eq. (B3) we derive the final expression for the 
intrinsic radiative rate of T-polarized QW excitons in the 
presence of homogeneous broadening: 



Tt{p\\ 



(V2k 2 p^) To 



(&q — p 2 ) 2 + 47 2 fc 2 p| 
1 



1/2 



(k 2 - p 2 ) 2 + Aj 2 k 2 P 2 - {k 2 - p 2 ) 



1/2 



1/2 



(B6) 
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For k -p\\ > 7/c Eq. (B4) is identical to Eq. (§g). The 
homogeneous linewidth Thom changes Eq. ( ftO|) for Tt(p\\) 
only at the band of states given by -y/co > k — p\\ > 
and, in particular, Tt(p\\ = k ) = T /(2^f 1 / 2 ) rather 
than divergent. Note that the band of states is rather 
narrow because 7 ~ l/(r x - x ujt) ~ 10~ 3 . The ID van 
Hove singularity at p^ = ko is absent for the L-polarized 
QW excitons. 

Homogeneous broadening practically does not change 
the optical decay rate of a gas of QW excitons. The use 
of Eq. (B6) for calculation of T opt by Eq. ( |3l| ) shows that 
for Maxwell-Boltzmann distributed QW excitons one ob- 
tains only a small renormalization of Jt of Eq. (J32[), i.e., 
Jt — * Jt(1 — I 1 / 2 )- As a result, the optical decay rate 
of a classical gas of QW excitons changes from T opt to 
r o p*[l — (3/4)7 1 / 2 ]. This correction due to the homo- 
geneous linewidth is very small and can indeed be ne- 
glected. The corresponding correction for statistically 
degenerate QW excitons is even smaller because the low- 
energy excitons tend to populate the ground-state mode 
P|i = rather than vicinity of the energy states E ~ Ek ■ 
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